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Abstract 

We report on calculations of the total gravitational energy radiated 
in the head-on black hole collision, where we use the geometry of the 
Robinson- Trautman metrics. 

PACS numbers: 04.30.Db, 04.25.-g, 04.25.Dm, 04.70.Bw 

1 Introduction 

Accurate description of the dynamics of relativistic gravitating systems remains 
one of the challenging problems in general relativity. There is great interest in 
particular on the description of systems involving very compact objects, as for 
example black holes. It is expected that these kind of systems will be describable 
in terms of an appropriate isolated model; which in a relativistic theory of 
gravity is expressed by an asymptotically flat spacetime. In this description the 
gravitational radiated energy will be encoded in the behavior of the physical 
field at future null infinity. 

In numerical calculations these ideas have been used in attacking the prob- 
lem of the time symmetric asymptotically flat initial data corresponding to a 
two equal mass black hole system[|l|||^. It would be nice to be able to compare 
these calculations with some analytic solution of the field equation; but in trying 
to do this one is faced with the fact that there is up to now only one known 
vacuum asymptotically flat radiating family of solutions, the so called Robinson- 
Trautman metrics!]^. These spacetimes represent the dynamics of a single black 
hole which settles down, in the asymptotic futureQ, to the Schwarzschild ge- 
ometry. 

In the process of two black holes falling into each other in a head-on colli- 
sion, one can distinguish two natural eras: one in which the two black holes are 
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separated and falling, and the other in which they already form a single black 
hole which is decaying to a stationary state. Since the Robinson- Trautman met- 
rics represent a black hole, one could estimate the total gravitational radiation 
emitted after the formation of the black hole using these geometries. In doing 
this one needs to associate the data of the two black hole collision system to the 
Robinson- Trautman data. We do this by matching, at an arbitrary retarded 
time, the radiation content of these metrics with the radiation calculated from 
the quadruple formula with Newtonian dynamics at the moment of collapse. 

2 Power of the energy radiated from Newtonian 
dynamics + quadrupole formula 

Let us consider the Newtonian problem of two particles of mass too, which start 
at rest at an initial distance di. One can estimate the power of the energy 
carried away by gravitational waves when they are at a distance d assuming 
Newtonian dynamics and the quadrupole approximation. More specifically; the 
quadrupole radiation formula relates the power of the energy radiated with the 
time derivatives of the quadrupole: 



and we are using geometric units in which the Newton's gravitational constant 
and the speed of light have the unit value. Then, the total power radiated is 



PQ{d) = 




where Mq = 2 toq is the total mechanical mass and d the distance at observation. 

3 Robinson- Trautman metrics: the quadrupole 
approximation 

In a previous paper Q it was found the general asymptotic behavior of the 
Robinson- Trautman metrics in the asymptotic future. The leading behavior of 
these metrics is governed by the quadrupole structure of the spacetimc. 




where 
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More explicitly, the whole geometry of these metrics is determined by a scalar 
V{u, 9, (j)) , where u, can be considered a retarded time and {9, (p) the coordinates 
of a 2-sphere. According to Ref. for an axisymmetric quadrupole structure, 
the asymptotic behavior of this scalar, for u — s- cx) is 

V = 1 + Ae^ -^¥20 + 0{e^ ^), 

where A is a constant Moo is the asymptotic mass of the spacetime in the regime 
u —^ 00, and Y20 is a spherical harmonic function. This geometry undoubtedly 
represents the dynamics of a single black hole; which becomes evident when in 
the limit for u 00 , i.e., V=l, one recognizes the Schwarzschild metric. 

As we have described above, there are two natural eras in the process of the 
head-on collision: before and after the formation of the remaining black hole. 
Since the Robinson- Trautman geometries describe the dynamics of a single black 
hole, it is natural to try to study, with these geometries, the final stage of the 
process. 

The time derivative of the Bondi shear of the Robinson- Trautman sections 
of future null infinity has the following asymptotic behavior: 



V6Ae "~ 2Y2o + 0{>), 



where now 2Y20 is a spin weight 2 spherical harmonic function and 0{>) means 
higher order terms of the time exponential. With this information one can 
calculate the flux of gravitational Bondi radiation at future null infinity in first 
order, obtaining: 

To relate the Robinson- Trautman geometry with the problem of the head-on 
collision of two black holes, we identify the Bondi power at some retarded time 
with the power calculated from the quadrupole formula at the moment when 
the two black holes touch; that is when d = 2Mq . In this way we determine A 
appearing in the above formulae. 

With the Bondi flux one can calculate in first order the total energy radiated, 
namely: 



E= r PB{u)du = 6^2M2o^4„/M„ 

Juo 47r 4 



47r 4 4 
Let Mi denote de mass of the spacetime at the initial retarded time; then one 
has the relation 

Afoo = M, - E. 

Therefore, the ratio radiation over initial mass, can be expressed by 
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Identifying the Bondi with the quadrupole flux one would obtain 




with If = ^ and I = . In this expression we identify M with Mi ; that is we 
map the Newtonian total mass to the total mass of the spacetime at the initial 
retarded time. 



4 Time symmetric 2-black hole initial data 

In reference it was presented the initial data corresponding to a time sym- 
metric 2-black hole system. There the ADM mass Mm was calculated, and a 
magnitude Dm (there denoted with capital 1) was defined, having the inter- 
pretation of some measure of initial distance between the black holes. More 
precisely Dm was defined as the length of the minimal line from the "ring" at 
minimum radius from one throat to the respective "ring" in the other throat. 

If one wants to compare calculations coming from different realizations of the 
same physical situation, namely the head-on collision of two equal mass black 
holes, one must provide a relation between the physical quantities involved. The 
two fundamental physical quantities involved in this situation are the total mass 
of the system and the separation of the black holes. As we have done in the 
previous section, the identification of the masses is straight forward, since the 
natural thing to do is to identify the invariantly defined total masses in each 
case. There remains however the relation among the notion of distance involved 
in each case. 

The question is: Is there a natural way to identify the Misner magnitude Dm 
with the Newtonian notion of distance appearing in the estimates of the total 
gravitational energy radiated using the quadrupole formula approach? A simple 
identification of Dm with d would be too naive in the regime of small distances. 
We seek then for some relation between Dm and d such that it captures the 
most important physical aspects of both systems. 

In order to do this let us recall that in reference the time symmetric initial 
data for a two equal mass black hole system was represented by the 3-geometry: 

where ds'^ is the line element of a flat 3-geometry, ri and r2 are the Euclidean 
distances from the field point to the deleted points in the fiat 3-geometry which 
represent the location of the first and second black hole respectively, and miao 
and m2oo are parameters which characterize the spacetime. If in the conformal, 
scalar appearing in the above line element, one eliminates the third term, the 
constant time Schwarzschild line element is obtained, with mass mioo expressed 
in terms of isotropic coordinates. 
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According to the arguments of reference Q the mass of, let us say, the first 
black hole is 

mi = mioo H 



2ri2 

where now ri2 is the Euclidean distance between the black holes; while the total 
mass of the spacetime is 

, , , mioc m2oo 
M = mi + m2 . 

This expression agrees with the first relativistic correction to the Newtonian 
total mass; since in first order one has 

2 

Mrei = Mo + Eo=mo + mo-'-^. 

di 

Therefore it seems that the Newtonian distance, for large d, should be identified 
with the coordinate distance of the flat 3-geometry of the time symmetric prob- 
lem. So using the line element ds^ shown above for the region between the black 
holes, and identifying the radius of each black hole with double their respective 
masses, we define the distance between the black holes by 

d— 2m2 , 
2m 1 

which in terms of the total mass M for a two equal mass system can be expressed 

by 
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d-A./(l + ^) 

then, defining L = ^ and I = jj one has the relation 

'-(1+3^) 



It is observed that with respect of this variable, the black holes do not touch 
exactly at I = 2 but when I* = 2 (l + jj^^ ; which corresponds to L — I* — 2.22. 
This picture is in agreement with the numerical calculations since it is found 
in reference Q that for L < 2.49 a common apparent horizon surrounds both 
holes Q. 

We relate d to the Misner magnitude Dm by setting D{d) = Dm- 

^The threshold mentioned in their reference is ^ < 1.36 which corresponds to our L < 2.49. 
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Figure 1: The total gravitational energy radiated away in the process in units 
of initial total mass. The two sets of numerical data correspond to the location 
of the detectors at L = 15 and L = 35 respectively. The vertical bar at approx- 
imately L — 2.49 corresponds to the region where a single black hole is formed. 
The RT estimate is plotted in terms of the relation L(l) explained in the main 
text. 



5 Upper limit for the radiation from the area 
theorem 

One can estimate an upper limit for the total energy radiated away in the form 
of gravitational waves when two very separated black holes collide, from the 
area theorem; which using the relation between total mass and single black hole 
mass shown above and in terms of the quantities already defined, establishes 
that 



6 Comparison of estimates with supercomputer 
results 

In ref. |^ the numerical calculation of the total energy radiated for the time 
symmetric head-on collision of two equal mass black holes was reported. They 
improved previous calculations 0| with the help of latest generation of super- 
computers, and further developed analytic and numerical techniques. In the 
next graph we plot their recent data, our Robinson- Trautman (RT) estimate 
for the total energy radiated and the upper limit from the area theorem. 
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Figure 2: Same as Fig. 1 but now replacing the plain RT estimate by a line that 
shows the total gravitational energy radiated away in the whole process, which 
is the sum of the RT energy estimate after the formation of the remaining black 
hole plus the quadrupole Newtonian estimate of the energy radiated before the 
formation of the remanent black hole. 



It can be observed that the RT estimate keeps bellow the supercomputer 
data; but this is no surprise since one had expected to only estimate, by this 
mean, the total energy radiated after the formation of the resulting black hole. 

The energy radiated in the first stage, when the black holes fall from di to 
d, before the formation of the resulting black hole, can be estimated from the 
quadrupole formula and using the Newtonian dynamics; that is 
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Eq = f ' Pq {x{t)) dt ^ I PQ{x)^dx- 



from which it is deduced that 

M~\l^^)\j\lf l) ^51 [if l) ^3P\lf I 

with I = j^, If = ^ and M is the total mass of the system. 

Then adding to the RT estimate the quadrupole Newtonian estimate one 
obtains the total energy radiated during the entire process in the quadrupole 
approximation, which is shown in the next graph. 

This remarkable agreement between our estimates and the supercomputer 
calculations allows one to think that the global quantities, as the total en- 
ergy radiated, for a gravitating system like the one under consideration, is not 
very dependent on the model realization of the system; since, in particular, the 
Robinson- Trautman approach is not time symmetric. Therefore the supercom- 
puter results obtained from the time symmetric two black holes head-on collision 
system might have a wider range of validity than the one assumed up to now. 
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It is important to note that our estimates can be calculated for values of L 
corresponding to distance over mass relation above the threshold of the forma- 
tion of the remaining black hole; which from reference |^ it is deduced to be 
around L=2.49. So our estimates are complementary of those calculated in ref- 
erence Q in which, using a completely different approach, they have estimated 
with high precision the total energy radiated for small L; while their method 
overestimates the total energy in two orders of magnitude in the vicinity of 
L fa 8 (which correspond to their /xq ~ 3 ). 

The total gravitational energy radiated in the large separation regime was 
represented through semianalytic method in reference |^ by their equation (6). 
The approach followed in that reference requires of a series of extrapolation in 
the applicability of known equations. First of all they use an equation for the 
total energy radiated calculated in the test particle limit|^ (eq. (1)); that is, 
one mass much smaller that the other mass ( m << M ) falling from infinity, 
and then they extrapolate the equation to the equal mass case {rh = M ). This 
approach needs from several correcting factors. Using an altered quadrupole 
radiation formula they suggest a correcting factor, appearing in their eq. (4), 
that takes into account the effects associated to the fact that the infall is not 
from infinity but from a finite distance. The alteration in the quadrupole formula 
comes from a substitution of the radial velocity by a "nonunique" |^ expression. 
It is observed that the integral in their eq. (4) is calculated up to the upper 
limit 2M , which is consistent with the test particle limit approach, but it 
overestimates the total gravitational energy radiation in the equal mass case. 
Therefore they have to introduce another correcting factor smaller than one, 
there called Ft, which is intended to take into account the black hole structure 
of both objects. These points were also studied in references and 
It is important to emphasize that we use, without alterations, the quadrupole 
radiation formula for the first stage before the collision, up to the point in which 
the two horizons touch each other ( d w AM ); this in contrast to the case of 
a test particle crossing one horizon ( d « 2M ). Then, in order to take into 
account the black hole structure of the remaining object we use the Robinson- 
Trautman geometry, which radiation is to be added to the quadrupole estimate; 
this in contrast to the reducing factor Fh . 

The present estimate of the total gravitational energy radiated is based on a 
couple of simple ideas with direct physical interpretation; which is a contribution 
to the understanding of the global aspects of the head-on black hole collision. 

It would be interesting to reproduce these calculations taken into account 
post-Newtonian techniques, as those used in reference , to see whether these 
effects will make any substantial difference to our calculations. 

A question that motivated us for this work was how a global quantity, as 
the total energy radiated, could be dependent on the details of the modeling of 
a physical situation. On purpose we use a completely different approach, and 
the result is that this global quantity does not appear to be very dependent on 
the details of the modeling. 

^See reference [14] of reference 
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As a final comment let us point out that since the time symmetric initial 
data obviously incorporates incoming radiation in the system, it is important 
to be able to know whether this incoming radiation will have an important con- 
tribution to the outgoing radiation. Since the Robinson- Trautman geometry 
can be considered to represent pure outgoing radiating spacetimes, and also 
the quadrupole estimate only takes into account retarded fields, the agreement 
of our estimates with the numerical calculations, of the time symmetric prob- 
lem, suggests that the incoming radiation contribution to the outgoing one is 
negligible. 
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